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We describe a proposal to probe the quantum tunneling mechanism of an 

individual ion trapped in a double-well electromagnetic potential. The time- 
evolution of the probability of fluorescence measurement of the electronic 
ground state is employed to characterize the single-particle tunneling mecha- 
nism. The proposed scheme can be used to implement quantum information 
devices. 
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Together with nonlocality and wave-vector collapse - the phenomena behind the recent 
advances in quantum communication [1] and computation [2] - the tunneling mechanism is 
one of the most intriguing aspects of the microscopic world of particles and their interactions, 
setting quantum reality apart from classical physics. While nonlocality and wave-vector 
collapse have been pursued as clues to understanding fundamental quantum phenomena 
such as the uncertainty principle and the process of quantum measurement [3], quantum 
tunneling is a valuable mechanism for probing the transition from quantum to classical 
dynamics [4]. 

Quantum tunneling at the macroscopic level has attracted much attention in the litera- 
ture [4-6], firstly owing to its application to SQUIDs (superconducting quantum interference 
devices), permitted by the advances in cryogenics, and recently because of the realization 
of Bose-Einstein condensation in dilute atomic gases. The observation of matter-wave in- 
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terference fringes has demonstrated that a Bose condensate consists of "laser-like" atoms 
which are spatially coherent and show long-range correlations, opening the field of coherent 
atomic beams and of atomic Josephson effect [7] . 

Probed at the single-particle level, quantum tunneling is just as provocative to our intu- 
ition as its manifestation on the macroscopic scale. The experimental techniques developed 
over the last decade for manipulating electronic and vibrational states of trapped ions can 
be employed to investigate fundamental quantum effects at a level so far accessible only as 
collective processes [8]. We describe in this letter a scheme for probing the tunneling mech- 
anism of an individual ion trapped in a double-well electromagnetic potential. The time 
evolution of the probability of a fluorescence measurement of the electronic ground state is 
used to probe the single-particle tunneling mechanism. 

We first assume i) that the ion is trapped in a symmetric double-well single-particle 
potential V(x) = b(x 2 — x1) 2 . (In ionic traps, the typical oscillation frequencies in the 
y and z directions are very much larger than that in the x direction, providing a good 
approximation to a one-dimensional trap.) The minima are given by Xq = ±^d/2b, so 
that the motion can be described as approximately harmonic, with frequency uj = ^Ad/m, 
m being the ionic mass. The harmonic approximation can be adjusted by fixing d (and 
consequently u ) and conveniently choosing the parameter b, which is equivalent to varying 
the height of the barrier separating the two wells, h = d 2 /4b. We also assume ii) that 
the parabolic approximation to the potential around each minimum is designed to contain 
(at least) the two lowest states of the harmonic oscillator described by the wave functions 



harmonic oscillator and i — 1,2 referring to both harmonic wells centered at x — =p£o- 
Finally, we assume in) that the potential is such that state <f>^ (of harmonic trap 1) is 
spatially close to state 0^ (of harmonic trap 2). Thus, these local modes are not orthogonal, 
due to the overlap e between the corresponding modes of the two wells: 




referring to the ground and first-excited states of the 




nm- 



(1) 
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When e< 1, given a first order correction , these local modes are approximately orthogonal 
and the eigenstates of the global double-well potential may be approximated by the sym- 



metric and asymmetric superpositions 4>±{x) ~ (x + xq) ± <f) 2 (x — x ) /V%, with 
corresponding eigenvalues E±^ = E^f 1 ± TZ^ n \ The coupling constants between the local 
states are given by the matrix elements 



n {n) = / dx 



^(x + x )\ [V(x)-V 1 [x + x ]\(f) 2 n >(x-x Q ), (2) 

where Vj [x — (— l) 1 ^] indicate the parabolic potential approximations around x = (— l) 1 ^. 
The tunneling frequencies between the corresponding eigenstates in V\ and V 2 are given 
by il^ = 2 \lZ^\ /h. In order to make the first order correction (e 1 ) valid, ensuring that 
TZ^/E^ < 1, we further impose the condition Xq/A^ 3> 1. Scaling the length in units of 
the position uncertainty in a harmonic oscillator ground state A 2 = fij (2muio), we obtain 
from Eq. (2) Q« = (3u x 2 / 8A 2 X ) exp (-x 2 /2A 2 x ) and tl^ = (x 2 /A 2 x ) 

The Hamiltonian describing the motional degrees of freedom of the ion trapped in the 
above- described double-well potential is 



motional 



r - h 2 

dx&H^, H = -—V 2 + V(x), (3) 
J 2m 



where \1/ (x) = 2~2 n i 4^ ( x ) * s the operator which annihilates the ion in the state 
4>[ n \ i.e., in energy level E^ of the harmonic well VJ. From the eigenvalue equation H(f)±^ = 
E^4^± an d the assumption that e< 1, the Hamiltonian in Eq. (3) becomes 

H moti onal = [Ejp (t^ + #H n) ) + ^ + t^)] ■ (4) 

Next, we assume that the trapped ion has two effective electronic states, excited |f) 
and ground ||), separated by frequency uj and coupled by the interaction with an effective 
laser plane wave propagating in the x direction, with wave vector = uj l /c. In this 
configuration, only the ionic motion along the x axis will be modified. The effective pumping 
laser beam, is detuned by 5 = uo—uol from the ||) <-> ||) transition. In NIST experiments with 
9 Be + the laser beams (with different frequencies) composing the effective pumping laser, are 



detuned from a third more excited level which, in the stimulated Raman-type configuration, 
is adiabatically eliminated [8]. In the Innsbruck experiments with 40 Ca + , a coherently driving 
direct transition between the electronic ground and excited states is employed [9]. A probe 
laser, strongly coupled to the transition between the electronic ground state and a third 
more excited level \s), is also employed in order to measure the ionic vibrational state by 
collecting the resonance fluorescence signal, which is the probability of the ion being found 
in the internal state ||). 

We can trail the tunneling phenomenon through the time-evolution of the probability of 
fluorescence measurement of the electronic ground state collected in two different situations: 

a) by individual addressing the ion in harmonic well 1 or, what seems to be more attractive 
experimentally, b) considering the laser beams (the pumping and the probe lasers) to reach 
both local potential wells simultaneously. The ion-laser interaction Hamiltonian that de- 
scribes the effective interaction of the quantized motion of the ionic center-of-mass coupled 
to its electronic degrees of freedom is [8] 

H ion - laser = % (?+ e^"*"^ + ~_ g-tffcx+iu^-i^ ? (5) 

where a + = |t) (II an d er_ = ||) (t| are the usual Pauli pseudo-spin operators, g is the 
effective Rabi frequency of the transition ||) <-> \\) and <fi is the phase difference between 
the two lasers composing the effective beam. Since, considering the situation a), the laser 
beams are disposed in harmonic well 1, the position operator x is given by 

x = y [ dx^f (x) x^ (x) cS n)t cS m) 

< *n,m J 

= A x J2 n (ci n)t c< n+1) + c^'cP) . (6) 

Evidently, considering the situation b) where the laser beams reach both local wells, the 
position operator turn to be x = A x X^=i 2 Vn + 1 (cf^ 1 cf" +V> + c^ +V> j . 

In what follows we consider only the situation a) since the generalization for the situation 

b) , which will be analyzed further, is straightforward. In a frame rotating at the effective 
laser frequency u>l, the ion-laser Hamiltonian is given by 
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Hi 0n -la S er = % exp (if] V U + 1 (c^V/^ + C^ 1 ^ C^) + ^i) + H.C.] , (7) 

where we have introduced the Lamb-Dicke parameter r\ = k^l\ x . 

In what follows we rewrite the total Hamiltonian in the interaction picture (script 
font labels) via the unitary transformation U(t) = exp ^—iH(,tj, where H = 
Yln-E^ +^"' ) ^2^) + indicates the free Hamiltonian composed of the in- 

ternal and motional degrees of freedom of the trapped ion, a z being the Pauli operator. The 
resulting Hamiltonian, including the motional and electronic degrees of freedom, reads 

n = J2 n K {n) fi n) % n) + &H n) ) + % exp (-t? 2 J2 n 3 B) M B) ) 

x {a + exp [^^^v/^^V^^e-^'^+^+^V^e^*)] e m+i ^ + H.c.} (8) 

In order to simplify the expression (8) we i) adjust the ion-laser detuning to the first red 
sideband for the ion-laser interaction (5 = u ). In addition, ii) the assumption of the 
standard Lamb-Dicke limit (77 <C 1), where the ionic center of mass is strongly localized 
with respect to the laser wavelength, enables the expansion of the Hamiltonian (8) to first 
order correction (77 1 ). Finally, Hi) with the optical rotating wave approximation, we obtain 
the expression 

n = Y, n nin) (% n H n) +#H n) ) 

+th V gexp (-rj 2 ^JC^) (*+ J2 n V^+T^f e +l ^ - H.c.) , (9) 

which in the absence of the tunneling mechanism leads to the coupling between the electronic 
and motional degrees of freedom of the trapped ion, described by the Jaynes-Cummings 
Hamiltonian (JCH). 

At this point, we assume that the ion is initially cooled to its motional ground state 
in harmonic well 1, and laser-excited to the electronic state Thus, the JCH induces 
the transition |1,0) 1 |f) <-> |0, l) l where the ket |1,0), (|0, indicates the ion in the 
motional ground (excited) state of V\. Evidently, being in the motional state |1, 0) 1 (|0, 1)^ 
the ion has a finite probability of tunneling to the corresponding state of V^: |1, 0) 2 (|0, 1) 2 ). 



Therefore, the motional basis states of the process are restricted to the ground and first 
excited states of the harmonic traps, which is why we assumed from the beginning that 
the parabolic approximation to the potential around each minimum is designed to contain 
(at least) the two lowest states of the harmonic oscillator. By collecting the resonance 
fluorescence signal in Vi, we can probe the tunneling mechanism through the behavior of 
the function Pj.(0 : the probability of fluorescence measurement of the electronic ground 
state. 

Considering the Lamb-Dicke limit and the fact that we have only one particle in the 
process, the function exp ^— ffYln^i^ ^1^) can be fairly approximated by unity, so the 
Hamiltonian governing the process has the simplified form (for n — 1, 2) 

n = (c?^ + c ( 1 1 ¥ 2 1)t ) + rc< 2 > (c^cf + <f<f t ) 

+ihrjg ( ( 7+c ( 1 1)t c ( 1 2) e +l ^ - (T-C^cfV^) . (10) 

Before proceeding further, it is worth mentioning that instead of adjusting the ion-laser 
detuning to the first red sideband for the ion-laser interaction (5 = u ), we could have chosen 
the carrier interaction where 5 = 0. In this case the effective Hamiltonian following from 
Eq. (8) (with the above-mentioned approximations) simplifies to 

1 ^carrier — /v I u l 4 T °1 °2 I ' /v \ 1 2 ' u l u 2 j 

+hg {a + e +llpL + o_e~ iLpL ) . (11) 

We note that the Pauli operators act only when the ion is in harmonic well 1. The dynamics 
governed by Hamiltonians (10) and (11) make it possible to probe the tunneling mechanism 
of the ion by collecting resonance fluorescence signals, as discussed above. 

Starting from the initial state |1, 0) 1 |0, 0) 2 ||) and Hamiltonian (10) we obtain from the 
Schrodinger evolution the evolved state 

hMt)> = [ci 1) (t)|i,o) 1 |o,o) 2 + c«(t)|o,o) 1 |i,o) 2 l |T) 



+ 



C[ 2) (t) |0, 1), |0, 0) 2 + Cf (t) |0, 0), |0, 1) 2 ||) , (12) 



whose coefficients satisfy the set of coupled linear equations 

ifrjC { i ] (t) = IZ^df (t) - (-l) n ihr)ge-^ ni ^C[ m) (t) , (13a) 
lUjCf } it) = U {n) ct ] it) , (13b) 

with n,m = 1,2 (n ^ m) and the initial conditions C[ 1] (0) = 1, c£ ] (0) = Cj 2) (0) = 

(2) 

C\ (0) = 0. Now, we introduce a further simplification into the present scheme. The larger 
the ratio Xg/A^ obtained by engineering the trap, the smaller becomes the ratio of the 
matrix elements \lZ^ /1Z^\. Therefore, when \lZ^ /1Z^\ C 1 we can neglect the tunneling 
process between the ground states in the two harmonic wells. In this case, the first term 
in Hamiltonian (10) can be disregarded and, by adjusting the phase of the laser pulse such 
that ifL = — 7r/2, we obtain the coefficients: 

C? (t) = [cos(£wt) + e - 1] /e, (0) = 0, (14a) 
C{ 2) (t) = i sin(£wt)/£, Cf (0) = [cos(£wt) - 1] /£ 2 , (14b) 

where the effective Rabi frequency, w = rjg, is modified by the parameter £ = 
< 1 + [TZ^ / (kqg)] > . Next we analyze the influence of parameter £ on the time-evolution 
of the probability of measuring fluorescence of the electronic ground state P^ (t) = 
harmonic well 1. Evidently, when VS 2 ^ = 0, we recover the well-known dynamics of -Pj.(t) 
for JCH [8]. For the choice TZ^/(hrjg) — 1 (£ = y/2) and employing the typical values in 
experiments with 40 Ca+ [9,10] r] « 0.1, g « 200 kHz, and oo « 2 MHz ( giving the estimate 
Tl {2) /(hr]g) ~ 150 x 1Z {2) /E^ ] ), we obtain 1Z {2) / « 7 x 10~ 3 as required to justify the 
approximation that the superpositions J (x) constitute eigenstates of the global double- well 
potential. In fact, remembering that TZ^/(hrjg) = (3uqXq/ lQr/gAl) exp (— Xq/2A 2 ), we ob- 
tain from the previous parameters the values x\j A 2 rs 17.3 and |7?. ( ' 1 - ) /7?. ( ' 2 ^ | ~ 6 x 10 -2 which 
are in agreement with the approximations considered above. In Fig. 1 we display the be- 
havior of function Pj,(t) for TZ^ — (£ = 1, corresponding to the dynamics of the JCH) and 
1Z^ I (hrjg) = 1. As anticipated by observing the expression for Pj.(t), the increase in the tun- 
neling rate IZ^ 1 leads to an increase in the effective frequency £w of population inversion and, 
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C { ? (t) 



m 



conversely, a decrease in the amplitude of the oscillations of Pj.(t), clearly indicating the tun- 
neling process. In fact, as soon as the ion reaches the excited state of local well 1, the coupling 
1Z^ to local well 2 entangles the motional excited states of both wells, preventing the proba- 
bility Pj.(£) (associated with the measurement of state |0, 1) 1 ||)) from reaching unity. When 
the curve in Fig. 1 for TZ^/(hrjg) = 1 reaches its maxima (wt = n rcy/2/ 4, n = 1, 2, ...), we 
have the entangled state \ |1, 0) 1 10, 0) 2 IT) + (75 |0, 1) 1 10, 0) 2 - \ |0, 0) 1 |0, 1) 2 ) ||). There- 
fore, both characteristics, the effective frequency £w of population inversion and the ampli- 
tude of the oscillations of Pj.(£), can be used to probe the single-particle tunneling mechanism. 

It is worth observing that the choice of the alternative initial state |0, 1) 1 |0, 0) 2 ||) leads 
to the entanglement (|1,0) 1 10,0) 2 |f) + |0, 0) x |0, 1) 2 |J,)) /y/2. In this situation, as expected 
for a closed system, we can re-establish a value of unity for Pj.(0- 

Now, turning to the Hamiltonian (11) and solving the Schrodinger equation for the initial 
state |0, 1) : |0,0) 2 we obtain the result 

p (f) = (Ai sin (Aigg - A 2 sin (X 2 gt)) 2 
lU (A 2 . -A 2 ) 2 ' 

where the parameters modifying the Rabi frequency, g, are Aj = 

r 1 1/2 1 1/2 

^ I 1 + 2 [R.^/(hg)} 2 + (-1) 1 1 + 4 [TZ^/(hg)} 2 \ . Evidently, when TZ^ = 0, we 
get the usual dynamics for the carrier pulse, with Pj,(i) = |sin {gt)\ 2 . In Fig. 2 we display the 
time evolution of Eq. (15), for TZ^ = and TZ^/(hg) = 1 (A; = [(3 + (-1)^) /2] V2 ). 
For the typical values given above, we obtain TZ^/E^ pa 7 x 10^ 2 (a value which match 
xl/Al 10.8 and \lZ^ /1Z^\ ps 0.1). Similarly to the behavior displayed in Fig. 1, we 
observe that the frequency of population inversion is higher than the Rabi frequency g as- 
sociated with free carrier dynamics. However, differently from the situation in Fig.l, the 
probability Pj_(t) (displaying the characteristic beat pattern due to Eq. (15)) can still reach 
unity. Since in the carrier regime the electronic states do not couple to the motional states, 
the excited and ground electronic states are both subject to the tunneling process (differently 
from the JCH case), resulting in the interference pattern shown in Fig. 2. 

From the value for the position uncertainty in a harmonic oscillator ground state A 2 , = 



h/ (2mw ) computed for 40 Ca + with motional frequency cj ~ 2 MHz, and considering the 
Jaynes-Cummings regime, where Xq/A^ 17.3, we obtain (using the typical values r\ 0.1 
and (7 ~ 200 kHz [9,10]) the distance between the local minima 2xo ~ 0.16 /im. Making the 
same estimates for the Carrier regime, where x^/Al ss 10.3, we obtain the distance between 
the local minima 2x « 0.13 /im. Recently, it was reported in Ref. [10] the laser addressing 
of individual ions in a linear ion trap with frequency 125 kHz when the distance between the 
ions is about 19. /im. The authors of Ref. [10] argue that the addressing technique presented 
permits individual addressing when the distance between the ions is only 7.6. /im with small 
error. Therefore, with nowadays technology it is difficult to individually address each local 
potential well. On the other hand, when considering the situation b) where the laser beams 
reach the local wells simultaneously, the Hamiltonian 10 (in Jaynes-Cummings regime and 
assuming \TZ W /1Z {2) \ < 1) turns to be 



f )e-^ . (16) 



Starting from the initial state |1,0) 1 |0,0) 2 ||) and Hamiltonian (16) we obtain from the 
Schrodinger evolution the evolved state described by Eq. (12) with the coefficients satisfying 
the coupled linear equations 

ihjC^ = ihge^CM (t) , (17a) 
ihjCW It) = K^C™ (t) - ihge-^cW (t) , (17b) 

with n,m = 1,2 (n ^ m) and the initial conditions C[ 1] (0) = 1, (0) = Cf ] (0) = 
(0) = 0. Solving the system above we obtain the coefficients: 

C? (t) = £ [X 2 ~ (A,) 2 + 3] cos (X j9 t) , (18a) 

CP it) = A* [(X 2 + 2) 2 - (X 2 + 2) (A,) 2 ] sin (X 3 gt) , (18b) 



p — i<-pL 



(x 2 + 4) j: l 2 



C ? ) W=^ Z ^E[ 2 ^) 2 -X 2 -2]co Sl 



(A^t), (18d) 



( r -i / >. 1/2 

where A , = 1 1 + x 2 + (-^Y X\/^ + X 2 /2} and x = ^ (2) / hrjg. The probability of 
measuring the electronic ground state ||) of the ion in both local wells 1 and 2 is given 

2 

by P^t) = C{ 2) {t)+C {2) (t) . In Fig. 3 we display the behavior of function P^t) for 
= (corresponding to the well know Jaynes-Cummings dynamics) and TZ^ / (hrjg) = 1 
( Ai = \J 3 — \fhj 2 and A 2 = \j 3 + y/E/ 2). Similarly to the behavior displayed in Fig. 1, 
we observe that the frequency of population inversion in Fig. 3 (full line) is higher than 
the Rabi frequency g associated with the free JCH dynamics whereas the amplitude of the 
oscillations of Pj.(0 is smaller than that for the free JCH dynamics. Both characteristics 
clearly indicate the tunneling process in Fig. 3, without the requirement of the individual 
addressing of each local potential well. 

In conclusion, we have presented a scheme for probing the tunneling process of a single ion 
trapped in a double-well electromagnetic potential. The tunneling dynamics is characterized 
by the behavior of the probability of collecting an electronic ground state fluorescence signal 
in one of the local potential wells. Two situations were analyzed, when collecting fluorescence 
measurement a) by individual addressing the ion in harmonic well 1, and b) considering the 
laser beams to reach the local wells simultaneously. Although the situation a) is beyond 
nowadays technology, the single-particle quantum tunneling can clearly be probed through 
the situation b). Considering the possibility of individual addressing of local wells, the 
present proposal can be extended to the implementation of the fundamental controlled-NOT 
two-bit gate by adding another laser beam in harmonic well 2 and storing the quantum bits 
in the motional states of both local wells [11]. We also note that our proposal can be 
implemented in spite of the decoherence processes coming from the coupling of the motional 
modes with the residual background gas [12] and with classical stochastic electric field [13], 
in addition to the finite lifetime of the electronic levels [13,14]. From the experimental results 
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reported in [8,9] we observe that several Rabi oscillations are visible within the decoherence 
time making it possible to observe the signature of the tunneling process. However, we stress 
that the designing of the double-well potential would require that the trap electrodes be only 
few micron away from the ion; then it is likely that the heating rate in the trap becomes 
higher than that observe in single traps, increasing the decoherence rate [15]. Besides, the 
required trap structure makes difficult to obtain high trap frequencies, about MHz, and 
consequently, makes difficult the cooling process of the ion to its motional ground state [15]. 
Here we note that the cooling process could be implemented in a single trap structure which 
could be adiabatically modified to a double-well potential. 

Although the design of such a double-well trap may turn out to be a considerable technical 
challenge, the fundamental principles discussed here can be implemented by engineering two- 
mode interactions in ion traps [16], in which the motional states of the ion in two different 
directions can be coupled, analogously to the dynamics for the double-well. Besides, the 
proposal here presented might provide a motivation for future experimental work. 
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Figure Captions 

Figure 1. The evolution of function Pj,(t), in the Jaynes-Cumming regime and individual 
addressing of local well 1, for the values TZ^/(hrjg) = (dotted line) and 1 (full line). 

Figure 2. The evolution of function P[(t), in the carrier regime and individual addressing 
of local well 1, for the values TZ^ / \hg) = (dotted line) and 1 (full line). 

Figure 3. The evolution of function Pj.(t), in the Jaynes-Cumming regime and the laser 
beams reaching both local wells simultaneously, for the values lZ^/(hrjg) = (dotted line) 
and 1 (full line). 



13 



